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Abstract 

In this letter we reconsider the previously given description of heavy mesons within a 
bosonized extended NJL model that combines heavy quark and chiral symmetry. In that 
work the naive gradient expansion of the quark determinant was used, which satisfactorily 
works in the light sector but does not adequately describe the heavy (0+,!"'") mesons. By 
investigating the exact momentum dependence of the quark loop we demonstrate that the 
naive gradient expansion in the heavy sector is not the right method to treat the unphysical 
gg-thresholds which would be absent in confining theories. We propose a modified gradient 
expansion which adequately extrapolates from the low-momentum region beyond threshold. 
This expansion gives a satisfactory description even of the (0^",1+) heavy mesons whose 
masses are significantly above threshold. 
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1 Introduction 



The observation of new symmetries for infinitely heavy quarks in QCD together with the for- 
mulation of the Heavy Quark Effective Theory (HQET) has led to a number of theoretical 
and phenomenological investigations. One interesting application is to exploit these symmetries 
in the formulation of tractable quark models ^, ^. 

In ref. Q we have combined chiral symmetry of light quarks with heavy quark spin and flavor 
symmetry in an extension of the Nambu-Jona-Lasinio (NJL) model. In order to be as close as 
possible to former investigations in the light sector (5|, |^ , we have fixed most of the parameters 
from light meson physics. For evaluating the quark determinant we have employed a linear 
gradient expansion. 

While our numerical investigations, relating heavy meson properties and model parameters, have 
been successful for the (0~,1~) meson doublet of heavy quark spin symmetry, our approach 
seems to fail for the heavier (0"^, 1"*") states of opposite parity. This seems somewhat surprising 
since such a problem does not arise in the analogous case of the p or ai mesons in the light 
sector. 

The aim of this work is to clarify this problem by a re-investigation of the gradient expansion in 
both the heavy and the light meson sector. The gradient expansion (or equivalently expansion 
of quark loops in external momenta) is usually performed, in order to avoid unphysical quark 
thresholds in NJL-type models without quark confinement. In fact, assuming such models valid 
in the low-momentum region far below the thresholds, a gradient expansion can afterwards be 
interpolated to external momenta of even a few 100 MeV above the threshold where it still 
predicts satisfactory results |5|, ^. 

Nevertheless, the situation for heavy and light mesons is qualitatively different due to the differ- 
ent analytical structure of quark loop expressions involving a heavy quark propagator {v-k+ie)~'^ . 
Therefore, one should look carefully how to define the interpolation procedure to larger external 
momenta which is necessary to describe the heavier states. 

The organization of the letter is as follows: In section 2 we will give a short review of the 
definition of the extended NJL model for heavy quark fiavors and of the result for the heavy 
meson self-energy which defines masses and renormalization factors of heavy mesons [Q]. We 
show that the (0~,1~) and (0'*',1'*') heavy mesons cannot be described simultaneously by one 
set of parameters. The next section 3 is devoted to a detailed analysis of the gradient expansion 
by evaluating the exact dependence on external momenta of the pertinent loop diagrams in 
both the light and heavy sector. Our analysis shows that a rather different way of interpolation 
is required for the heavy meson self-energy as compared to the light sector. This suggests a 
modified gradient expansion which is valid in both the light and heavy sector. The modified 
gradient expansion is shown to lead to the desired improvement in numerical results. Finally, 
some concluding remarks are given in section 4. 

2 Free effective meson lagrangian and lowest— order gradient ex- 
pansion 

In ref. |^ we have presented an extension of the NJL model which combines chiral symmetry 
for light quarks with heavy quark symmetries for heavy quark fields defined by 

Qv{x) = ^—^eyiY>{imQV x)Q{x) . (1) 

We do not want to give technical details of the employed bosonization procedure but simply 
quote the result for the effective meson lagrangian when all the quark fields have been integrated 



1 



k 



H,K 




H,K 



k — p 



Figure 1: Self-energy diagram for heavy meson fields H,K. 
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where A'^c = 3 is the color factor, rfiQ is the current mass matrix of light quark flavors and Gi, 
G2, G3 are coupling constants of four-quark interactions. Furthermore, 



iI^ = i^-T. + \f + 4-f^-{H + K) {iv ■ dy\H + K) 



(3) 



is the Dirac operator for the light quarks, which contains the several light (S, V^, A^) and heavy 
meson fields {H,K). In the heavy quark limit, heavy meson fields are organized in spin symmetry 
doublets 
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with $ being the heavy scalar, <I>^ the heavy pseudoscalar, the heavy vector and <^^^ the 
heavy axial vector field, respectively. V, A denote the light vector and axial-vector fields. For 
the light octet of Goldstone bosons we use the common non-linear representation ^ = exp^iir / F) 
where vr = ■K°'\'p/2 and F is the bare decay constant. Vector and axial-vector expressions are 
build via = i/2{id^e} + S^d^i), Al = i/2{id^i^ - (}d^i). Finally, the light scalar field S 
achieves a non-vanishing vacuum expectation value indicating the spontaneously breaking of 
chiral symmetry. 

From (H) we have derived self-energy expressions for mesons as well as interaction terms between 
heavy and light mesons. This is achieved by evaluating Tr Inii/* in terms of quark loops that are 
regularized by Schwinger's gauge invariant proper-time method. Note that due to chiral and 
heavy quark symmetry, the heavy sector is controlled by only one coupling constant G3 for both 
H and K fields, while all other parameters are fixed from the light sector. 

The self-energy term for heavy mesons ^h,k as a function of the external momentum {v ■ p) is 
represented in figure || and calculated as 



-tro [H'll^Hiv ■ P)H'\ + trn • p)K' 

reg ^4^ iYjiU^ - i) + m^) (H' +T) {H^ + 
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where i = u,d, s is a hght flavor index. 

Remember that the heavy (Q^g)-meson lagrangian in the heavy quark limit is given by^ 



^heavy ^ 

+trD 



H\iv ■ d - AMh)H' 

K\iv d - AMk)K' , (7) 



where AM|^^ = M\j — mq, with M|^^ being the heavy meson mass. Let us compare the 
inverse propagator term (v-p — AMjj j^) resulting from (^) with the self-energy expression 
of (^) and the term in dH). First observe that I^^h x ^ function oiv-p only. Expanding 

^H,K around the on-shell value v ■ p = AMj^ one obtains 

nk,/f (t^ • P) = irH^K{^MH,K) + IIIk{AMh,k) {v-p- AMh,k) + 0{{v ■ p - AMh^k?) . (8) 
With this expansion one reads off from (|6|) that the meson mass is determined by 

li{AM\j^j,) {AMh,K ± m*) + /j - ^ = (9) 

and the Z-factors for the necessary field renormalization H^,K^ {Z\j W ^ K'^ are given 

by 

Z'H,K{^Mh,K) = {ll{^Mh,K) + '^li{^Mh,K) (AMI,,^ ± m')y' , (10) 
where !{, are loop integrals defined by 

J 1 - 



167r4 J P - {m^y 
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^(m^)2r(-l,(m*)VA') , (11) 




(A;^ — (m*)^)(f ■ k + V ■ p + ie) 







^ {m^y-{v-pY 



^^{m^Y-[v.pfT\-\ r > I i> , (12) 



I^iv ■ p) = 1/2 dl^{v ■ p)/d{v ■ p) and r(Q, x) is the incomplete gamma function. 
As a consequence of the heavy flavor symmetry the solutions of (^) do not depend on the heavy 
quark mass. But they do depend on the light quark constituent masses and so are influenced by 
the explicit and spontaneous chiral symmetry breaking. For that reason we have supplied the 
mass shifts AM\j ^ with the index i of the light quark flavors. 

It is obvious that for v-p > the quark loop produces an imaginary part due to the unphysical 
quark-antiquark threshold. Clearly, one must neglect this imaginary part of the meson self 
energy, since it would yield a nonvanishing decay width T(H^,K^ ~^ 'f + Q) forbidden by 
confinement. Its formal appearance is obviously related to the fact that confinement is not 
taken into account in our model explicitely. The simplest way to avoid these unphysical two- 
quark thresholds is to expand Uj^ ^ to first order around (w • p) = in analogy to the gradient 

^ Note that the heavy meson field has the mass dimension 3/2. 
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AMfj - AM^''^ [MeV] 


240 


80 


20 


<0 


100 


Md., - Md [MeV] 


AM^"^ - AM^'^ [MeV] 


2710 


1450 


1030 


820 


410 


Md, - Md [MeV] 


AM|. - AM]^'" [MeV] 


9020 


4780 


3370 


2670 


530 


Md^, - Mr) [MeV] 



Table 1: Differences of the mass shifts as a function of the heavy-Ught quark couphng constant 
G3. The column "Exp." contains the experimental values. 



expansion of the light meson sector, where one expands to first order around = [||, ^]. This 
lowest-order expansion in (v ■ p) has been performed in ref. [^], leading to 

AM|,,^(/^(0)±2m*/^(0))±m^I^(0) + /i-^ = , (13) 

ZH,K = {im±2m'li{Q)y^ . (14) 

Note that eq. (^) is linear in AM, and the Z-factors in eq. (14) are independent of AM. 
Let us also recall the numerical results following from our simple gradient expansion ( [l^Jl^ ). 
While the values of the several mass differences AM\j ^ are to be compared with experimental 
values, the Z]:^ ^-factors enter into the decay constants fH,K of heavy mesons defined by 

(0|g7^(l-75)Q„|//40~)) = ifnMHV^ , 

(0|g7/.(l-75)Q.|i^t,(0+)) = -IkMkv^ , (15) 

via the equation Q 



fH,K./M^=^^^^ . (16) 
The light parameters m''''^ = 300 MeV, = 510 MeV, A = 1.25 GeV, entering (|l|,|l|) are 



fixed from light meson properties, whereas the heavy-light quark coupling constant G3 has to 
be adjusted from heavy meson physics. 

As table |l] shows, it is impossible to find a consistent value of G3 that fits both, the masses 
of heavy H and K mesons simultaneously. While a value of G3 = 8.7 GeV^^ gives realistic 
results for masses and decay constants of the members of the H{0^,1~) multiplet [^], the 
masses and Z-factors for the parity conjugate mesons K come out simply too large. This is 
in contrast to the light fiavor sector where both vector and axial vector mesons are adequately 
reproduced. Obviously something must go wrong in the gradient expansion in the heavy sector. 
We therefore abandon the naive gradient expansion and re-investigate the exact expressions 
(^^). This will lead us to a modified gradient expansion which overcomes the shortcomings of 
the naive gradient expansion in the heavy sector. It is important to note that such a modification 
of the gradient expansion should not be understood as a better approximation (recall that we 
even know the exact result) but rather as a heuristic method to extract physical information 
from quark loop calculations without including unphysical threshold effects. We hope that the 
following discussion will illuminate these ideas. 



3 Modified gradient expansion 

Let us now start to investigate the exact equations (PJIOj). Since we expect the mass shifts 
AMj^ to be larger than m* we have to handle the imaginary part of the integral (|l^ ) above 
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the threshold. For illustration let us drop the imaginary part and show the real part of eq. (P). 
In figure 2 we have plotted, for convenience, Gg ^ as a function of AM]^ resulting from the 



exact solution (^) and from the simple gradient expansion (13), respectively. (Heavy mesons 
with strangeness show a completely similar behaviour.) Note that already for values below 
the treshold AM = m the linear approximation seems problematic. Instead, G3 ^ as a function 
of AM seems to start already with a significant curvature. Keeping more of the quadratic 
behaviour at AM = would immediately decrease AMk for a given value of G3 considerably, 
whereas AMh would not such dramatically change. 

For clarification let us compare the renormalization factors Zjj- ^ in the exact formula ( p!o|) with 
the simple gradient expansion (p!^). As an example we plotted {Z'fj)~^ in figure 3. Once again 
the simple gradient expansion cannot be viewed as a good low-energy approximation. This 
becomes even more appearent if we consider the analogous situation for the light sector. In 
figure 4 we plotted the real part of an exact calculation of the renormalization factor for the p 
meson given by 

z-p\v') = Y^l dxix{i - x)r(o, ^2 — —) ■ (17) 

Here the usual gradient expansion is indeed justified as a definite interpolation from the low 
energy region to regions above the threshold. 

Let us now try to understand analytically which terms are responsible for the deviation from 
the linear behaviour of eq. @ (respectively from the constant behaviour of the Z-factors) in the 
low-energy region. The threshold presents itself in the combination m? — (x)AM^ as one sees 



from the integral Isdv ■ p) = AM) in eq. (12). For AM > m the incomplete gamma function 
and the square root with negative argument will produce imaginary parts. In the proper-time 
formalism this is the usual way how thresholds arise. Since the low-energy region is defined as 
AM <C m we are on the safe side if we neglect AM^ compared to m? in the integral expressions, 
performing the replacement 



Is (AM) ^ -^{2 AMY 




Nc 

= /3(0) + 2AM/2(0) , (18) 
J2(AM) ^ /2(0) . (19) 



Indeed such a strategy would reproduce a constant Z-factor for the p meson. Furthermore, for 
heavy mesons the integral is multiplied by (AMitm), and this will generate a quadratic term 
(in AM) in the self-energy. For light mesons such an additional pre-factor is missing. Note that 
our so defined modified gradient expansion has not to be viewed as an ordinary Taylor expansion 
of the self-energy. The crucial point of this investigation is that it would mean no improvement 
to calculate any higher order terms in the gradient expansion since we would get more and 
more unphysical effects from the threshold. Instead, our method represents a compromise for 
extracting relevant physics from the low energy region and neglecting unphysical imaginary parts 
and threshold effects leading to the emission of free quarks. 

One could try to treat the thresholds in the quark loop integrals seriously and simply dropping 
imaginary parts. For the light (axial) vector mesons this would (accidently) make not much 
difference. However for heavy mesons such a strategy gets no numerical support and seems 
unphysical at all. 

With this understanding we can dare to extrapolate from the region below threshold to mass 
shifts AM above threshold. Doing so we hope that chiral symmetry dominates the physics over 
a large region, at least so far that also the (0+, 1+) multiplet can be described by our modified 
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Figure 2: The behaviour of I/G3 as a function of AM]^/m" (bold line) and of AMf^lm'^ (thin 
line). The relation (13) representing the lowest order gradient expansion is plotted in each case 
with the same line type and can be identified as a linear function originating from AM = 0. 
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Figure 3: The behaviour of the normalized (Zfj)~^ as a function of AM]^/m" (bold line). The 
threshold at AM/m = 1 shows up as a singularity. In lowest order gradient expansion the 
Z-factor stays constant at the value at AM = (thin line). 
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Figure 4: The inverse normalized Z-factor (Zp 



as a function of M^/(4m^) 



(bold line). This 



time the threshold at 



4m^ shows up as a peak. The lowest order gradient expansion (thin 
line) is for a large region in good agreement with the exact solution. 
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gradient expansion. Equations (|9D, (|10D are then to be replaced by 

1 



2Gs 



(20) 
(21) 



Figure 5 shows the behaviour of ^3"^ in the modified gradient expansion (20) to be compared 
with figure 2. 

In table |2| we present the numerical results following from this modified gradient expansion. 
Note that heavy quark masses are extracted from m^^c = Mb,d — AM]j. For the calculation 



G3[GeV-^] 


3 


5 


7 


9 


11 


"Exp." 




nic [MeV] 


1070 


1410 


1600 


1740 


1840 


1550 


Mj,^/2 [MeV] 


mb [MeV] 


4470 


4810 


5000 


5140 


5240 


4730 


Mt /2 [MeV] 


AM^ - AM^" [MeV] 


250 


160 


110 


80 


50 


100 


Md, - Md [MeV] 


AM]^''* - AM^'" [MeV] 


370 


390 


410 


430 


450 


410 


Md^ - Md [MeV] 


AM|, - AM]^''' [MeV] 


870 


800 


780 


780 


780 


530 


Md^, - Md [MeV] 


fB [MeV] 


300 


210 


170 


150 


130 


180 


1 


>i 


fu/fH 


1.13 


1.13 


1.12 


1.11 


1.09 


1.H.2 


i 


,1 



Table 2: Here we show different observables as a function of the coupling constant G3 calculated 
with the modified gradient expansion. The column "Exp." contains the experimental values or 
in case of the heavy constituent quark masses some plausible estimates. In the last column we 
give the way of determination. (For the electroweak decay constants we have given the values 
calculated from the lattice and QCD sum rules.) 

of heavy meson decay constants we have used the experimental values of meson masses if 
available, or have made use of the relation Mb^ — Mb = Md^ — Md which is true in the heavy 
quark limit (x = {s, 1, si}). 

Indeed we observe a drastic improvement compared to the results in table For a value G3 ~ 7 
GeV~^ a simultaneous fit of nearly all heavy meson properties is possible. Only for the most 
massive {Bsi,Dsi) states, masses are still predicted slightly too high^. 

With our analysis we are now able to predict the decay constant of several B-mesons (Note that 
fD is assumed to get large l/rric corrections). Within a physical acceptable range of 

6.5 GeV-2 < G3 < 7.5 GeV'^ (22) 

we obtain 

160 GeV < fs < 180 GeV (23) 

180 GeV < fB, < 200 GeV (24) 

150 GeV < fs, < 160 GeV (25) 

155 GeV < < 175 GeV . (26) 

^ This indicates that a simple local interaction is perhaps not sufficient to describe a heavy quark weakly 
bound via a p-wave with a not so light strange quark. 



9 




0.5 1 1.5 2 2.5 3 

AM/m 



Figure 5: The behaviour of I/G3 as a function of AMj^/m" (bold line) and of AMf^lm'^ (thin 
line) using the modified gradient expansion (20) and the exact solution (9), respectively. 
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4 Conclusions 



In this paper we have studied in detail the treatment of quark loop diagrams, arising from the 
calculation of the quark determinant in the recently proposed extended NJL model 
Usually, the problems originating from the lack of confinement in quark models which leads to 
unphysical quark thresholds, are overcome by an expansion in external momenta. The lowest- 
order gradient expansion works successfully for composite light mesons. However, in the heavy 
meson sector there is phenomenological evidence that the gradient expansion has to be improved, 
which has been discussed in section 2. 

We have shown that this fact is due to an essentially different analytical structure of the loop 
integrals under concern, coming from the special form of the heavy quark propagator {v-k+ie)~^ . 
While the renormalization factor Zp for p mesons stays nearly constant below the threshold, this 
is not the case for the corresponding expressions Zh,k of heavy mesons. We have convinced 
ourselves that this additional momentum dependence is not due to threshold effects and should 
be included. We stress again that our modified procedure should not be viewed as an ordinary 
Taylor expansion that would never be able to interpolate beyond the threshold. 
Most importantly, with this improved gradient expansion heavy meson properties of both, the 
(0^, 1^) and the (0^, 1^) spin symmetry doublet can now be described simultaneously with a 
quark-coupling constant G3 ~ 7 GeV~^. With this value we are in a position to estimate several 
weak decay constants of B mesons. 

In summary, it has been shown that the evaluation of quark loop diagrams with heavy quark 
propagators in unconfined models demands a modification of the gradient expansion that has 
been used in the light sector. A phenomenologically successful procedure seems to consist in 
neglecting unphysical two-quark threshold effects in a definite way and keeping all remaining mo- 
mentum dependence. Phenomenologically, this leads to realistic results for heavy meson masses 
and decay constants within our extended NJL model with chiral and heavy quark symmetries. 
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